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Abstract 



The two Higgs doublet model for the top (T2HDM) is a model with two scalar doublets in 
which the top quark receives a special status. The special status of the top is manifest in the 
Yukawa potential, by coupling it to the second Higgs doublet, while all other quarks couple to 
the first Higgs doublet. The working assumption of the model is that the vacuum expectation 
value (VEV) of the second Higgs (v 2 ) is much larger than the first Higgs VEV (i>i), so that the 
top receives a much larger mass than all other quarks in a natural manner, and tan/5 = i^/fi 
is large. In addition, these Yukawa couplings generate potentially enhanced flavor-changing 
(FC) interactions, both in the charged and the neutral sectors. These interactions can greatly 
enhance FC decays such as t — > ch and h — * tc. 

In this work we explicitly (and independently) derive the Yukawa and Higgs potential of the 
T2HDM, obtaining the scalar to quarks and triple scalar interactions Feynman rules. We 
calculate the branching ratio (BR) of the one-loop and tree-level rare FC decays t — > ch and 
h — > tc in the T2HDM. We explore the BR within the parameter space of the T2HDM, focusing 
on regions in which BR {t — > ch) and BR {h — > tc) in the T2HDM can be enhanced compared 
to these BR's in the standard model (SM) and two Higgs doublet models (2HDM) of types I 
and II. We find that the BR of the rare decays t — > ch and h ^>ic can be enhanced by many 
orders of magnitude in the T2HDM compared to the BR in the SM and in the 2HDM-I,II, 
especially in regions of the parameter space where the decays are dominated by dynamics of 
the neutral scalar sector. 

The BR (t — > ch) can be measured in the upcoming Large Hadron Collider (LHC), if its value is 
above ~ 5 x 1CT 5 which is the discovery threshold of the LHC. We find that the BR (t — > ch) can 
exceed the LHC threshold in certain regions of the parameter space of the T2HDM, reaching 
up to ~ 10~ 4 . Moreover, we find that the BR(h — > tc) in the T2HDM can exceed ~ 10~ 4 . 
Discovering these processes at the LHC will show a clear indication of new physics beyond the 
standard model, and will particularly motivate the special dynamics of the T2HDM setup. 



Chapter 1 
Introduction 



The standard model (SM) of elementary particles has been highly successful in describing 
observed and measured phenomena. It contains, however, an unexplored sector, namely, the 
Higgs sector. The SM also has several problems, one of which is the quark mass hierarchy 
problem, especially the top quark having a much larger mass than all other quarks. 
In its minimal form the SM Higgs sector is comprised of one Higgs doublet, but that is not 
necessarily the case. Non minimal extensions of the Higgs sector can describe the same observed 
phenomena, and predict additional phenomena, which are as yet unobserved, but not ruled out. 
This work will describe one such extension of the SM - the two-Higgs doublet model "for the 
top" (T2HDM). The T2HDM features particular Yukawa couplings whereby the top quark 
receives a special status. This particular Yukawa structure also gives rise to potentially large 
flavor-changing (FC) couplings in the up-quark sector. 

In this work we will explicitly (and independently) derive the Yukawa potential of the model, 
though it has been shown elsewhere [1, 2]. 

The FC rare decays t — > ch and h — > ic have a very low branching ratio (BR) in the SM, 
of ~ 10" 13 [3, 4]. This low BR makes these decays extremely sensitive to new physics in the 
scalar sector. In this work we will explore the BR of the FC rare decays t — > ch and h — > ic 
in the parameter space of the T2HDM, at the 1-loop level (we adhere in this work to the 
t'Hooft Feynman gauge) and at the tree-level order. We will focus on regions of the parameter 
space in which the BR (t — > ch) can exceed the detection limit of the upcoming large hadron 
collider (LHC), and also on regions where the BR(t — > ch) and BR(h — > ic) can be enhanced 
significantly compared to other two Higgs doublet models (2HDM). 

1.1 The two Higgs doublets model 

The minimal extension of the SM is the two Higgs doublet model (2HDM). A comprehensive 
review of the principles of the 2HDM can be found in [5]. Basically, the model is comprised of 
two Higgs doublets, <3>i and <3>2- They usually obey discrete symmetries, whose aim is to define 
the Yukawa terms, and which divide them into several types. The 3 most common types are: 
type I, where $ 2 couples to all quarks, and $i does not couple to quarks; type II, where $i 
couples to down quarks, and $2 couples to up quarks; type III, which denotes a general case 
in which both $1 and $2 couple to all quarks. 

The type II 2HDM describes the Yukawa structure of the minimal supersymmetric standard 
model (MSSM), and is therefore of particular interest in the literature. 

Several properties are common to all types of 2HDM. All feature additional physical scalars: 
initially there are two complex doublets, hence 8 (2x4) degrees of freedom, and the electroweak 
breaking absorbs 3. We are therefore left with 5 degrees of freedom which are equivalent to 5 
physical scalars, plus the 3 (unphysical) Goldstone bosons which are present also in the minimal 
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SM, and which are "eaten" by the gauge bosons. The components of the scalar doublets mix to 
produce the mass eigenstates. They are denoted as follows: h°, H° - CP-even neutral scalars, 
A - CP-odd neutral scalar, G° - neutral (unphysical) Goldstone boson, If ± - charged scalars, 
G^ - charged (unphysical) Goldstone bosons. The mixing conserves the symmetries of the 
theory: the mass matrix does not mix scalars with different charges, and in CP conserving 
theories there is no mixing between CP-even and CP-odd scalars. 



1.2 The two Higgs doublets model "for the top" 

The 2HDM "for the top" (T2HDM) was first introduced by Das and Kao [1] as an effective 
approach for providing the top its mass in a natural way. They proposed a 2HDM in which 
the second Higgs field couples only to the top, while the first Higgs field couples to all other 
quarks. 

The choice of the coupling can be expressed also in terms of a discrete symmetry imposed on 
the Lagrangian [1], under which the fields transform as follows: 



$i -> d R ^-d R , u R ^ -u R (u = u,c), 

$2 -> +$2, Ql -> +Q L , t R -> +t R (1.1) 

where: are the Higgs fields, (u, c, t) R are the right-handed SU(2) singlet up-type quarks, 
d R = (d, s, b) R are the down-type right-handed quarks, and Ql is the left-handed SU (2) quark 
doublet. The discrete symmetry (1.1) produces the Yukawa couplings of the T2HDM, described 
below. This discrete symmetry is softly broken by the A5 term of (2.2) as discussed in [6]. 
With these Yukawa couplings, the top gets its mass primarily from the second Higgs vacuum 
expectation value (VEV), which we will choose to be much larger than the first Higgs VEV: 

->1. (1.2) 

Vi 

This is the working assumption of the T2HDM. 

This particular Yukawa coupling can also give rise to large FC interactions, as we shall later 
show. 

Distinct features of the T2HDM are: 

• The H + cb vertex is enhanced by the ratio of CKM matrix elements V t b/V c b compared to 
other 2HDM's. This property motivated our work, as well as the analysis in [7, 8]. 

• There are tree-level FC interactions in the up-quark sector; but there are no tree-level 
FC interactions in the down-quark sector, unlike the case of the 2HDM-III in which the 
tree-level FC interactions are both in the up and down-quark sectors. 

• The couplings of the neutral scalars (if , h°, A ) to all the quarks except for the top quark, 
increase with tan/3. This property motivated the analysis in [9]. 

These points will be further elaborated upon in Sec. 3. 

The T2HDM could stand on its own, although the couplings and symmetries defined above 
do not seem "naturally derived". However, it could also be viewed as an effective low energy 
realization of a more fundamental theory. Some examples are: 

• An extra-dimensions scenario, Randall-Sundrum like, in which the couplings are derived 
from the location of fields in the 5th dimension [10]. 
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• A technicolor scenario with a topcolor condensate scalar having a large VEV, which by 
construction couples only to the top quark [11]. 

• A non-minimal supersymmetry scenario, in which <3>i couples to down quarks and $2 to 
up quarks, but the couplings for ur, Cr are very small, and get most of their value from 
loop corrections. 

1.3 Rare processes and the t — > ch and h — > ic rare decays 

Rare decays are a sensitive probe for new physics [12]. Such decays are defined as rare because 
in the SM they are subject to a suppression mechanism, which can be either highly effective 
cancellations, such as the GIM mechanism, or the conservation of a fundamental symmetry, 
such as lepton number. 

New physics models may greatly enhance such processes, by working around the suppression 
mechanism. For example, the process \i — > cy which is forbidden by lepton flavor conservation, 
can be realized by relaxing the symmetry in the neutrino sector (see e.g. [13]). 
In this work, we have chosen to explore the BR of the rare decays h — > ic and t — > ch in the 
T2HDM, as a potential new-physics signal at the LHC. The LHC discovery limit for the t — * ch 
decay process is BR > 5.8 ■ 10~ 5 [14] for an integrated luminosity of 100/6 -1 . As mentioned 
before, the SM BR is about ~ 10~ 13 , and, therefore, unobservable at the LHC. Previous studies 
[15] have shown that the BR (t — > ch), where h = H°, h°, A , could reach up to ~ 10~ 4 in the 
2HDM type II and in the MSSM, and about ~ 1(T 6 in the 2HDM type I. 
In the T2HDM the FC decays can be enhanced due to the large H + cb coupling which is 
proportional to V t b x tan (3 instead of V c b x tan (3 in other 2HDM's, as we shall later show. This 
large coupling motivated us in calculating the BR of t — > ch. 

As an aside, we will briefly recall how the experimental detection of the process will proceed 
at the LHC [14]. At the LHC the top will be mainly produced in ti pairs. One then searches 
for processes in which the t decays to ch, while the i decays in the main bW~ channel. The h 
is most likely to decay into bb pairs when its mass is below 130 GeV, whereas above this mass 
the W + W~ decay channel starts to dominate [16]. The full process (in the lower mass range) 
will look like: gg — > ti — > hcbW~ — > bbcblv. The main background will come from a similar 
process in which the t — > bW + — > bjc (where j denotes a quark jet), and the i decays as before. 
In this case a misidentification of the jet as b will result in an erroneous t — > ch identification 
[14]. 

The h — > tc decay is the complementary process to t — » ch if mh > m t + m c . The amplitude 
of the process is equal to the amplitude of t — > ch, by applying crossing symmetry [17], and 
therefore it is subject to the same enhancements as the t — > ch process, compared with other 
2HDM's and the SM. 

1.4 Predictions and constraints on the T2HDM 

To date several rare decays and other observables have been calculated in the T2HDM: 

• The electric dipole moment (EDM) of the electron was calculated in the T2HDM [1], and 
the neutron EDM in [18], for their dependence on the CP violating mixing in the Higgs 
sector. The experimental results constrain this mixing. 

• The process b — > was calculated in the T2HDM in the leading order, for its contribution 
to 67,8 [19, 20]. By adding this result to the SM prediction one can compare the theory 
to the experimental result: BR{b — > sj) = (3.55 ± 0.26) x 10~ 4 [21], and derive bounds 
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on the model. A prediction was also given for the partial rate asymmetry of the decay, 
which is very different from the SM prediction, and can be measured in B factories. Newer 
measurements in [21] seem to further restrict the additional CP violating phase in the 
Yukawa sector of the T2HDM (see below). 

• The meson mixings B — B [22], K — K [8], D — D [8, 1] were calculated for contributions 
to the mass splittings ArriB,K,D, £k of K — K mixing, the ratio p/q of D — D mixing. 
These results further constrain the parameter space of the T2HDM, in a manner similar 
to b — > S7, as was discussed above, ruling out different regions of the parameter space 
[1,8,22]. 

• The process b — > sl^lj was calculated in [23], and was found to constrain the T2HDM 
weakly, so that those bounds are included within other calculated bounds. 

• The process gq — > qqq was considered in [9], where q denotes a b or c quark. It was found 
that the T2HDM cross section for this process will be detectable at the LHC, while the 
MSSM and the 2HDM type II are not expected to have a detectable signal. Therefore if 
such a signal is observed at the LHC then it will stand out as a clear indication in favor 
of the T2HDM. 

• The process Z — > bs + bs was calculated in various models [7]. Experimentally the BR 
has a weak upper bound. It was found that the T2HDM BR for this decay is comparable 
in size to the SM predicted value (~ 1CT 8 ), and to the MSSM with t — c mixing. In 
comparison, MSSM with b — s mixing is about two orders of magnitude higher, while the 
2HDM type II is about two orders of magnitude lower. 

• In a recent article [2] some of the above calculations were simultaneously combined for a fit 
to recent experimental data, mostly from B-factories. The best-fit values and la intervals 
for all the parameters in the fit were calculated. As this is the most comprehensive work 
constraining the T2HDM parameters, these were the bounds used in the present work. 

The rare decays t — > ch and h —> ic have not been calculated yet in the T2HDM. This work 
is aimed at this calculation, with the intention of giving a prediction which will hopefully be 
verifiable at the LHC. 
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Chapter 2 

Yukawa interactions in the T2HDM 



In this section we give an explicit derivation of the Feynman rules of scalar-quark-quark inter- 
actions in the T2HDM. We will start from the interaction-basis Lagrangian, which follows from 
the symmetries imposed. We will rotate the quark fields and the scalar fields to their mass 
basis. Finally, we will write the Yukawa terms in the mass basis, arranged by interactions, in 
terms of standard parameters. 

The Lagrangian density of the T2HDM Yukawa interactions is of the following form [1]: 

C Y = -QiA&iFijdRj - Qu^iG ij=1 ,2 ^ " ^ - QiA^iG a tR + h.c. , (2.1) 

where: i,j = 1,2,3 are flavour indices, L(R) = (1 — (+)7 5 ) /2 are the chiral left (right) pro- 
jection operators, fi{R) = L(R)f are left(right)-handed fermion fields, F,G are general 3x3 
Yukawa matrices, and: 




The Higgs potential can be generically written as (assuming CP conservation) [5]: 

C H = Ai ($+$! - vl/2) 2 + A 2 ($+$ 2 - vl/2) 2 + A 3 - v\/2) + ($£$ 2 - v 2 2 /2)} 2 + 

+A 4 [(^$i) K$ 2 ) - (^$2) + A 5 \$t*2 ~ viv 2 /2\ 2 . (2.2) 

The absence of CP violation implies that the CP-even and CP-odd Higgs mass-eigenstates do 
not mix, and that the VEV's can be taken to be real without affecting the Lagrangian of the 
theory [5]. 

Dropping the flavor indices and defining: I^ 12 ) = diag(l, 1,0), 1^ = diag(0, 0, 1), the Yukawa 
potential reads: 

C Y = -Q L $iFd R - Q L foGlM + $ 2 GI< 3 >) u R + h.c. . (2.3) 
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Inserting $, $ into the Yukawa potential and rearranging the Yukawa terms: 
C Y = -d L ^=Fd R - u Li ^= ( Vl GlW + v 2 GI^) u R - Q L ^Fd R - 

-Q L (^GI^ + $ 2 GlW) u R - ^Qu$iG l3 t R + ^Q L AG i3 t R + h.c. = 
= -d L ^=Fd R - u L -j= faGI™ + v 2 Gl^) u R - Q L ^Fd R - 

—Ql ($iGI (12) + ^iCI {3) ^) u R - Q Lt (® 2 - ^$1) G l3 t R + h.c. = 
= -d L ^=Fd R - u L ^= ( Vl GI^ + v 2 GI^) u R - d L -j=Q\Fd R - u L ^Fd R + 

_ % _L$0* G ^(12) + ^j(3)^ Ur + ^-q ^(12) + ^j(3)^ Ur + 

- Ul ^= ($°* - ^$;*) ci< 3 >f * + d L ($ 2 - - oi (3) ^ + h.c. . (2.4) 

If the Gij are of 0(1), and v 2 is much larger than vi, then the eigenvalues of the matrix 
(viGI^ + v 2 GI^^j can be expanded as a series of V\jv 2 . After expanding to the leading order, 
the mass matrix eigenvalues are: 0(1) • [vi, v±, v 2 ]. As can be seen, the top quark receives a mass 
contribution from the second, and larger, VEV, while the up and charm receive their masses 
from the first VEV. 

Rotating to the quark mass basis, we define: di,R — > F>L,RdL,R, ul,r — > Ul,rUl,r, such that: 



M d = ^j=D\FD R = diag (m d , m s , m b ) , 

M u = U{^= (vrGl^ + v 2 Gt 3) ) U R = diag (m u , m c , m t ) . 
V2 



(2.5) 



We define the CKM matrix: Vqkm = U^Dl, Vq KM = D^Ul (we will henceforward drop the 
subscript CKM when referring to the CKM matrix), and a new mixing matrix for the up-quarks: 

E ee M U U R I^U R , (2.6) 

as was originally defined in [1]. 

The matrix U R can be generally parametrized by multiplying 3 rotation matrices [8]. Defining 
3 rotation angles: a i2 = 0, a 23 = sin" 1 (e ct £) and a± 3 = sin" 1 (e ct £'), where e ct = and 

are parameters naturally of 0(1), and £ = \£,\ tiPi , we get: 



/ cos — sin 
sin cos 
\ 1 



/ 1 



/ 



\ / V 1 - l^f o -e^'* \ 

^/l-le^l 2 -egg 1 _0 

(2.7) 



\ 



where the asterisks (*) denote terms which are not relevant for our calculations to follow. 
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Using Eq. (2.7) we can now write the E matrix: 

I "^A l^'l 2 i 1 ~ \^) ^4^1 - Ml* (l - \e ct tf) yfl - | 2 

V ^' (i - M 2 ) ^/i-kctei 2 eet^-I^V 1 -! 6 ^'! 2 ( x - l £c ^l 2 ) - l e ^'l 2 ) 

(2.8) 

The special mixing matrix between up and down quarks via the charged Higgs, neglecting terms 
of O {e 2 ct ) or O (^^j (recall: e ct = ^), can then be written approximately as: 

( e ct ev td ew*ty cs +e ct ev ts e ct ev tb \ 

(sV)/m f = a% el t \Z\ 2 V cs + e ct CV ts e ct CV tb . (2.9) 

V V* $£Vo + V ts V tb - V th el (|e| 2 + in 2 ) / 

We note that the matrix was approximated for the purpose of illustration. In all calculations 
the matrix was used without neglecting anything. 
The Yukawa terms in the quark mass basis are: 

C Y = -d L MdR-u L M u u R -h®\—d R -u L $¥—u R - 

Vi Vi 

-u L Ul^D L D[FD R d R + d L D[^U L U[- ( Vl Gl™ + v 2 GI^) U R u R - (2.10) 

Vi 

-u L Ul^= ($? - Gl^U R u R + d L D{ - U L UlGl^U R u R + h.c. . 

We will henceforward drop the mass terms. 

Using: I( 12 ) • IP) = [0] 3x3 , and: U[Gl^U R = U* L ± ( Vl GI^ + v 2 GI^) U R U R I^U R = ^E, we 
get: 

Cy = -d L <$>\^d R -U L <$> Q *^U R -<$>tu L V- -d R + ^d L Vi KM - ±U R - 

Vi Vi Vi Vi 

_ _ V A$>A u L -Hu R + ( $ 2 - - ^) d L V^Eu R + h.c. . (2.11) 

The Higgs fields are not in their mass basis yet. The mass basis is derived from the Higgs 
potential in the Lagrangian, see App. A (see also [5]). 

We note that the real parts of the neutral Higgs fields generally rotate with the angle a, whereas 
the other Higgs fields rotate with the angle f3. This is due to their shift by the VEV's. This 
point is sometimes overlooked in the literature, where often the value for a is chosen arbitrarily. 
This omission of a can perhaps be attributed to the MSSM, where the angle a is constrained 
by the Higgs masses [5]. 

The Higgs fields in the mass basis are defined as follows: 
$° = $ 0r + tan (/?) = — , 

$° r = if cos a sin a, $f = G° cos (3 - A Q sin ft, $+ = G + cos/3 - H + sin/5, 
$° r = # sina + /i°cosa, ^ = G° sin [3 + A cos (3, $+ = G + sin/3 + H + cos (3. 



(2.12) 
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The Yukawa terms in the quark and Higgs mass basis, are then: 



C Y = -d L [(H°cosa-h°sina) +i(G° cos (3- A sin (3)] ^±d R - 
-u L [ (H° cos a - h° sin a) - i (G° cos (3 - A sin 0) ] ^u R - 

- (G + cos /3 - if + sin /3) u L V^^d R + (G-cos/3-H-fwP) d L V^^^u R - 

- [ (H° sin a + 0° r cos a) - i (G° sin (3 + A cos 0) ] u L —Zu R + 



v 2 



+ — \(H°cosa-<f) 2 r *sma) - i (G° cos (3 - A sin (3)] u L —Zu R + 

fl V 2 



+ 



(G~ sin (3 + H cos 0) - ^ (G - cos /3 - if " sin /3) 



d L ^ f — Zu R + h.c. . (2.13) 

t>2 



Collecting identical interactions: 



C Y = -d L [(H° cos a - h° sin a) + i(G° cos (3 -A sin p)]^d R + 



+H°u L 
+h°u L 
+tG°u L 
+iA°u L 



M u L 1 „ 

cos a E sin a H E cos a 

Vi v 2 v 1 



1 



1 



sin a E cos a E sin a 



v 2 



Vl 



cos p H E sin p E cos p 

v 2 Vi 



u R + 
u R + 



— - sin (3 H E cos /3 H E sin (3 

Vi v 2 Vi 



u R 



(-G+ cos (3 + H + sin 0) u L V^^d R + 

v ' Vi 



+G~d L V* 
+H-d L V^ 



cos p h sin p E — cos p E 

Vi v 2 Vi 



u R + 



— sin p h cos p E + sin p E 

Vi v 2 Vi 



u R + h.c. . (2.14) 



Using: v\ = vcos[3, v 2 — vsin(3, v — \Jv \ + v | = 2i!1hi j anc i adding the h.c, we get the final 
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Yukawa terms, in the physical mass basis, arranged by interactions: 



C Y =H"d 

+ A°d 

+ H°u 

+ h°u 

+ A°u 
+ G + u 



+ H+ 



u 



gM d cos a 
2mw cos (3 

~.gM d 
i- — 

2m w 
9 

2m w 
9 

2m w 

'. 9 
i- 

2m w 

9 

\[2m w 
9 

\[2m w 



d + h°d 



tan (3{R-L) 



gM d sin a 

2mw cos (3 

d + G°d 



—i- 



d+ 

gM d 



. , cos a _ 
-M n + E 



cos (3 



2m w 
sin a cos a 



(R — L) 



d+ 



+ 



sin (3 cos (3 



R + (h.c.)L 



u+ 



„ sina _ 
M„ - E 



cos (3 



cos a sin a 
+ 



i?+ (/i.c.)L 



sin /3 cos (3 
(-M u tan/3 + E (tan (3 + cot (3)) R + (/i.e.) L 

[-yM d i? + M U FL] d + h.c.+ 



u+ 

u + G°u 



2m w v 1 



u+ 



[tan (3VM d R+ (-M u tan/3 + E t (tan/3 + cot /?)) VL] d + h.c. . (2.15) 



For completeness, in App. B we give the complete list of Feynman rules for the T2HDM model, 
as derived above. 
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Chapter 3 

The flavor-changing sector of the T2HDM 



The T2HDM features unique flavor-changing (FC) couplings of both charged Higgs and neutral 
Higgs, where the neutral Higgs FC couplings are in the up-quark sector only. These couplings 
can enhance FC processes, in particular the t — > ch decay. 



3.1 Charged-Higgs FC Yukawa interactions 

In the SM the leading order diagrams of the t — > ch decay are at the 1-loop level, through 
the mediation of W + gauge bosons and b-quarks in the loop, being oc V cb . In general, in any 
2HDM there is a corresponding vertex of H + cb, and the process t — > ch proceeds also via similar 
diagrams with H + and b in the loop. In the T2HDM the H + cb vertex receives a particular 
value distinct from other 2HDM's: 

L D — H + c [tan (3VM d R + (-M u tan (3 + E f (tan (3 + cot (3)) VL] cb b ~ 
v 2mw 

~ —J- — H + c [tan (3V cb m b R + m c (- tan (3V cb + f * (tan (3 + cot (3) V tb ) L] b . (3.1) 
\J1m w 

As can be seen, the H + cb vertex has a term proportional to V cb x tan/3 which is common to 
other 2HDM's, but has an additional term proportional to (tan/3 + cot (3) x {Y)V} cb ~ m c ^*V tb) 
as shown above in (2.9). The main contribution for 1-loop diagrams with internal H + and b, 
in the T2HDM, will therefore come from m c ^*V tb terms, and is thus not CKM suppressed. 
For our analysis we also consider the H + ib vertex in the T2HDM: 

L D -jM- — H + i {tan f3V tb m b R+ 
V2m w 



+ 
9 



m t 14tan/3 + m t (y tb - V tb e 2 ct (\£\ 2 + |ef)) (tan/3 + cot /?)] l} b = 
H + i {tan (3V tb m b R+ 



+ 



m t V tb cotp-m t V tb e 2 ct (\i\ 2 + |^'| 2 ) (tan (3 + cot l} b (3.2) 



We can see that by taking — > the H + tb interaction in (3.2) becomes equivalent to that 
of a 2HDM type I or II. 



3.2 Neutral-Higgs FC Yukawa interactions 

A priori there is no distinction between h° and H° other than the rotation angle a. In this 
work we therefore adopt a = (3, and in the following discussion we explore the consequences of 
this choice, in particular the elimination of the tree-level H°ic vertex. 
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The t — > cH° decay in the T2HDM can proceed at tree level, from the following interaction 
term (see Sec. 2): 



L D HH 



9 



2m w 
9 



, , cos a ^ 
-M u + E 

cos p 

sin a cos a 
+ 



2m w \ sin p 1 cos {3 



sin a cos a 
sin (3 cos p 1 

(m c ^i? + m c e ct iL) 



R+{h.c.)L 



c ~ 



tc 



(3.3) 



where we used the off-diagonal terms of E from Eq. 2.8, neglecting terms of order e 2 ct (recall 
e c t = m c /m t ), S te ~ m c £ and (E t ) tc ~ m c e ct £. 

For arbitrary a and f3, this will lead to t — > cif° decay at tree level. On the other hand, the 
H°tc vertex can vanish if: 
1)£ = 0, 



2 -- 



sin a i cos a 



sin /3 cos 8 



a = f3 + n7r . 



In this work we wish to examine the case of a vanishing H°tc tree-level interaction, therefore 
adopting a — (3, since: 

1. £ = is strongly disfavoured by the bounds in [2], as we will discuss in Sec. 4.2. 

2. £ = cancels the potentially enhanced term of the H + cb coupling, as shown above. 

3. The limit a = f3 is a natural result of the MSSM, when the mass of the CP-odd neutral 
Higgs, A , is large (see e.g. [5] in the limit m A o — > oo). Note that a = f3 is widely used 
in the literature, partly for this reason. Even though the T2HDM setup is not natural 
within the MSSM, this will help us compare our results with other existing results in 
different types of 2HDM's. 

4. The limit a = (3 sets the scalar H° to be SM-like. As such, it will have SM-Higgs Yukawa 
couplings: 



£y =/3 D H°d 



gM d 



2m 



w 



d + H°u 



gMu 
2m 



w 



u . 



(3.4) 



In that case, the direct mass bounds on the SM Higgs may roughly apply to H°. 



We turn now to the /i°-up quarks Yukawa interactions. The 1-loop t — > cH° can proceed also 
through the mediation of h° scalars and top quarks in the loop. In these diagrams the important 
interactions are the h°tc and the h°it vertices, which get special values in the T2HDM. 
The h°ic interaction reads: 



L D h°t 
~ h°t 



9 



2mw 
9 



. , sin a _ 
M u 

cosp 

cos a sin a 
+ 



cos a sin a . . 

+ ) ) R+(h.c.)L 



2m w \ sin (3 cos (3 
where setting a — j3 gives: 



sin (3 cos (3 

m c £R + m c e ct £L) 



c ~ 



tc 



C, 



(3.5) 



C a=(} D h°i 



9 



2m w 



(tan (3 + cot f3) (m c ^R + m c e ct ^L) 



c. 



(3.6) 



One can see that the h°tc interaction can be enhanced by tan/5. 
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The h°tt interaction, with a — /3, reads: 



L D h°t I — — [m t tan/3 - E tt (tan + cot 0)\ R + (h.c.) L \t ~ 

-m t cot /3 + m c e ct (j£| 2 + |£'| 2 ) (tan /3 + cot /?)] i? + (/i.e.) l| i, (3.7) 



where we use S tt = (S^) ~ m t — m c e ct (|£| 2 + |£'| 2 ) . As in the case of H + tb, the leading term 
oc m t tan/3 cancels, leaving terms that are suppressed either by e 2 a or cot 2 /5. 
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Chapter 4 



Calculations 



4.1 One-loop amplitude 

The 1-loop decay amplitude is composed of 10 Feynman diagrams, shown in Fig. 4.1. Their 
explicit calculation is given in App. C. 

The calculation was aimed to be model-independent. Thus, the Feynman diagrams in Fig. 4.1 

were drawn by assuming general vertices, which are defined in Fig. B.l in App. C, and by 

assuming general fields: 

5j- denotes a quark (up or down type) 

V a - denotes vector (gauge) fields 

H a denotes scalar fields 

In this way it was possible to calculate the same process in different models, by inserting the 
appropriate vertices and fields. 

The 1-loop integrals were calculated numerically with FORTRAN using the FF package [24]. 
The calculations were done using the Passarino-Veltman reduction scheme, which expresses the 
integrals in terms of basic scalar n-point functions. All other (vector, tensor) integrals can be 
computed using combinations of the scalar functions (for explicit furmulae see e.g. [15] App. 
A). In App. D we give the definitions of the reduced functions that we used to calculate the 
1-loop integrals. 

4.2 Bounds on the parameter space of the T2HDM 

As was mentioned in the introduction, bounds on the T2HDM charged-sector parameter space 
were simultaneously calculated in [2] to give a best fit to various experimental results. The 
processes that were selected were the ones most sensitive to the T2HDM. 
We describe here the processes for which the parameter fit was done in [2], and summarize in 
Eq. (4.1) the allowed values of the parameters of the T2HDM at la. 

• The BR of B — >• X s ^ [20] was estimated from the shifts in the Wilson coefficients Cr,s 
caused by the T2HDM. The BR constrains mainly the parameters m#+, tan/3 and £. 

• The BR of B + — > t + u t was calculated in the T2HDM [2]. This process receives a tree- 
level contribution from charged-Higgs exchange, and has a large impact in constraining 
the parameter space, especially m H + and 

• The CP- violating parameter Ek was calculated in the T2HDM [8], and was found in [2] 
to severely constrain the parameters m#+, tan/3 and £. 

• The time dependent amplitude of the CP asymmetry cl^k = A (_B° — > J /ipK^j, was found 
to constrain the parameters m H + and tan/3. 
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t — > ch diagrams 
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• The neutron EDM constrains mainly m^, tan/3 and 

• The Amp mass difference from D — D mixing can be completely dominated by NP 
effects, and so it was required in the analysis that the T2HDM value would not exceed 
the experimental value. This requirement constrains m H +, tan/3, £ and 

• The ratio Am^/Am^ was included in the fit, although it gave weaker constraints than 
e K . 

In the analysis of [2], a \ 2 function was defined, that featured as variables the T2HDM charged 
Higgs sector parameters - tan/3, mjj, f and £' and the SM CKM-matrix parameters - p, 
r) and [a,/3, 7] (unitarity triangle). These parameters were simultaneously fit to the processes 
described above. As this is the most comprehensive work constraining the T2HDM parameters, 
these were the bounds used in the present work. Since they are directly relevant to the present 
work, we list the final results of [2] below (recall that f = |£| e tlpi ): 



m H ± 


= (6601118) GeV, 


tan/3 


= 28lf, 


lei 


~ 0.8, 0.5 < |f | < 1 




= (no±gg)°, 


\a 


~ 0.21, 




~ 250°. 



We note that these values are also allowed for the 2HDM-II as given in [25]. We further note 
that the Amp mass difference from D — D mixing receives a contribution from the neutral 
sector that can constrain the neutral sector parameters, but that contribution is suppressed by 
/ m \ 2 

a factor of ( ^^^7"^- J compared to the charged Higgs contribution, and therefore does not 

impose significant constraints on the parameter space of the neutral sector [26, 8]. 
Direct constraints from experiments impose weak bounds [27]: for the SM Higgs, the direct 
search bound is m H o > 114 GeV. For supersymmetry, bounds for neutral scalars are m h > 
90 GeV, while for charged scalars current bounds are m H + > 80 GeV. 



4.3 From amplitude to BR 

From the amplitude we get the width of the decay using [27]: 



r = 4,.Ai(l,4,4).L, (4 . 2) 
where 5^|.M| 2 is the squared amplitude averaged over initial polarizations and summed over 

pol 

final polarizations, and A (x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. In the case of an incoming 
fermion, we have: E\M\ 2 = ^El-^l 2 - 

pol pol 

Calculating the squared amplitude, we get: 
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Ill 

M = T7^(M L L + M R R)u t , 

l07T 



J2\ M \ 2 = 7^— 4 tr[u c (M L L + M R R)u t u t (M* R L + M* L R)u c ]= (4.3) 



256tt 4 

pol 



256tt 4 
The BR is then: 



[2m c m t (M L M* + M R M* L ) + (m 2 c + m\ - m 2 h ) (M L M* L + M R M* R )] 



Where usually only x = W + b is taken, since its BR is very close to 1 [27]. We take the tree-level 
value of r (t — > W + b) = 1.55 GeV. If the mass of the H + is smaller than the mass of the top, 
the process t — > H + b is possible and significant, and must be taken into account in the sum. 



4.4 Higgs decay BR 

The 1-loop amplitude of the opposite process, h — * ic + ct, is identical to the 1-loop amplitude 
of t — > ch, by applying crossing symmetry [17]. The width of the process is given by: 



2 2X EM 

r (h - tc+ ct) = 2 x r (h -> *c) = 2iV c At ( 1, ^| ) • ^ , (4.5) 



where iV c = 3 is the color factor, and Xll-^| 2 = Yl \-M\ 2 in this case (an incoming scalar). 

pol pol 

In order to calculate the BR for h — > ic, one has to calculate the total width of the scalars. 
In this work we have included leading-order contributions to the Higgs width from h — > qq, 
h — > VV, h — > 2 scalars and /i — >■ vector+scalar [28]. The last can be important in some 
regions of the parameter space such as low tan (3. The formulae used for the calculation of the 
total width are given in App. E, along with a plot of the SM Higgs width in the leading order 
approximation compared with higher order predictions, in Fig. E.l. 



4.5 Checks of the calculations 

• We have successfully reproduced the results for BR(t — > ch) obtained in [3] in the SM 
and in [15] in the 2HDM-II. 

• We have successfully reproduced the results for BR (h — * ic) obtained in [4] in the SM 
and in [25] in the 2HDM-II. We note that Arhrib's results matched for a E w ~ 1/128.9, 
which is different from the one reportedly used, a E w ~ 1/137, as also confirmed by him 
in a private communication. However, we were not able to reproduce the BR(h — > ic) 
values of [15], as also stated in [25]. 

• Some amplitudes have a divergent part. Since the process is calculated at leading order, 
no renormalization is needed to cancel the divergent terms, and they should cancel among 
themselves. This cancellation is demonstrated in App. F analytically. It was also verified 
numerically in the FORTRAN code. Checking that the results do not diverge is also a 
test of the self-consistency of the 1-loop amplitude calculations. 
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4.6 Tree- level amplitude 



As stated in the introduction, there are FC tree-level interactions of h° in the T2HDM, when 
a = P, as opposed to H° which has no FC tree-level interactions when a — (5. Therefore, the 
decays t — > ch° and h° — > tc at tree-level are possible when allowed kinematically. The tree-level 
coupling h°ic was given in Sec. 3. Here we derive the leading order tree-level BR values of 
these two decays. We will neglect throughout the derivation terms of order: m 2 /m 2 , m 2 /mf, 
m c / m h an d cot 2 p. The last term is neglected in accordance with the working assumption of 
the T2HDM, which is a large tan/3. We also set a — p. The last requirement is not imperative, 
but it renders simpler formulae and makes the following derivation consistent with the 1-loop 
calculations of t — > cH° and H° — > tc decays in this work. 
The tree-level amplitude for the process t — > ch° is: 



M (t -> ch°) 



u, 



9 



2m w 

where from (2.8) we have: 



. , , . sin a 
(M u ) ct = - S 



cos/3 



cos a sin a , . 
a I —15 + — 15 ) )R + (h.c.)L 



sin p cos P 



(4.6) 



S rf = m t e 2 ct C\ll - kctCTV 1 ~ l e ct£T ~ m c£ctC 



(St) ct = - \e ct ^\ 2 y/l - \e ct ?\ 2 ~ m c t, 



(4.7) 



which we insert in the amplitude with a = /3 to get: 



M(t—> ch°) = u c —^- (cot P + tan/5) m c f [e ct i? + L] ct u t 
2my/ 

= u c [M R R + M L L] ut. 



(4.8) 



The squared amplitude summed over external spinors is: 



\Mf = 2m c m t (M L M* R + M R M* L ) + (m 2 t + m\ - m 2 h ) (M L M* L + M R M* R ) = 



g 2 m 2 
Am 2 w 



2 \y-\2 



(cot/? + tan/3) |f 



[2m c m t ■ 2e ct + (m 2 + m 2 c - m 2 h0 ) (l + e^) ] 



~ 9 ^t^p\t?[rn 2 -rn 2 hQ ]. 



(4.9) 



From Y,\M\* = |E \M\ 2 we can calculate the BR's of the processes t — > ch° and /i° — > tc. 

po/ po( 

The width of t — > c/i° reads: 



r (t — > c/i°) = 47r • f 1,4.^1 



g 2 \£\ 2 m t m 2 c 
\2%ixm 2 w 



tan 2 /3 1 - 



poi 



QAir 2 mt 



m 



h° 



mi 



(4.10) 



The width for t — ► 6VF + (at tree-level and neglecting terms of order m\jm 2 and a s ) is [27]: 



r (t -> w + ) 



3V 

64vr 



1 - 



m 2 



1-2 



m 2 



+ 



m 



w 



(4.11) 
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from which we get the leading order BR (t — > ch°) (for large tan/3): 

BR(t ch") ~ l -p£t^P (l - 4V fl - (l -24 + 4V- (4-12) 

2m^ V m t J \ m t J \ m t ™w) 

For instance, for the best-fit parameters of Eq. (4.1), tan/3 = 28, |£| = 0.8, and for a — (5, and 
= 91 GeV, we get BR(t -> c/i°) = 0.0077. 

1 1 2 

By applying crossing symmetry on , the tree-level squared amplitude of the opposite 
process, h° — > tc, is: 



,2^2 



El^l'-f^^^l^Ko-rn?], (4.13) 



from which we get: 



/ 2 n El-^l 

r - tc + ct) = 2 x r (h° - tc) = 2N C \^ ( 1, ^r, ^J- j poi 



^0 "^0 

7Vc| l l / 2 T m - tan^fl-4V- (4-14) 
327rm^ \ m h ° J 



lQirrrihO 



-',.2 • / ^2\2 





For a — P, (assuming also that m^o < 2m^o, 2m H +) the total width of /i is mainly comprised 
of fermion decays, since the couplings W + W" h° , Z°Z°h° and H°H h° are all oc sin (/3 — a) 
(see table B.2 and App. E). Below the tt threshold (at about 340 GeV) the bb decays dominate. 
The width of h° — > bb is then (see App. E): 

In this case, the BR of h° — > tc is: 

Bi?(/,^tc+ct)^|e| 2 ^(l-^) 2 . (4.16) 

m b V ''V/ 

For |f | = 0.8 and m h o = 300 GeV, we get BR{h° -> tc+ a) = 0.023 . 
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Chapter 5 
Results 



In this section we give our results and discussion for the 1-loop decays BR (t — > cH°) and 
BR (H° — > tc) in the T2HDM. All the masses are in units of GeV. We set a = f3 for the reasons 
explained above. Other parameters are set to their best-fit value of (4.1) unless stated otherwise. 
Our calculations were done in the t'Hooft Feynman gauge. The SM Higgs has a best-fit mass to 
EW precision data of 9ll 3 2 GeV [27]. In our setup in which H° has couplings identical to the 
SM Higgs, we expect these bounds to be roughly applicable, and therefore we set m H o = 91 GeV 
in the process t — > cH°. We take the total t-quark width r (t — > W + b) = 1.55 GeV. For the 
process H° — > tc we arbitrarily choose m H o = 300 GeV. 

For definiteness, other values we used for the calculations were [27]: m t = 172.5 GeV (pole 
mass), m c = m c (m c ) = 1.24 GeV, m b = m b (m h ) = 4.20 GeV (m c and m& are in the MS 
renormalization scheme), mw = 80.40 GeV, mz = 91.188 GeV, cos9w = rnw/mz, chew (jn z ) = 
1/127.9. The values were used as given here, without running in energy scale. The BR results 
were found to be sensitive to m c and m&: for example, for the BR value 5.99 x 10~ 5 quoted in 
table 5.1 in the upper row for the T2HDM with m b = 4.2 and m c = 1.25, setting m b {mz) ~ 3 
and m c {m z ) ~ 0.7 [29] yields 1.28 x 10~ 5 . We note that our results were not sensitive to m s . 
We used m^o = 1000 GeV to enhance the triple-scalar coupling which is roughly oc m 2 A0 , as can 
be seen in App. A and B. The £ and £' parameters are set to their best-fit value of Eq. (4.1): 
|f | = 0.8, ipz = 110°, |f | = 0.21, and (p? = 250°. 

5.1 Results for the 1-loop top rare decay t — > cH° 

In Fig. (5.1 a) we show a 3D plot of BR{t — > cH°) in the m H + — tan/3 plane in the T2HDM. 
The flat grid in Figs. 5.1 and 5.2 is the LHC detection limit of BR > 5 ■ 10~ 5 , so that the 
colored surface above the grid is the region in the parameter space which has (in the T2HDM) 
a BR high enough to be detected at the LHC. 

The choice m h o = 1000 GeV in Fig. 5.1 suppresses the diagrams in which the neutral h° Higgs 
runs in the loop and, thus, better explores the charged Higgs sector properties. As expected 
the BR rises with tan/? and is highest when 17Ih+ is lowest. The dominant Feynman diagram 
in this case is the one which has two H + scalars and a b quark in the loop, and is shown in Fig. 
(5.1 b). This diagram receives an enhancement from the 3-scalar vertex, as discussed above. 
In Fig. 5.2 we show the BR in the m^o — tan/5 plane in the T2HDM. 

We took m H + = 1000 GeV in Fig. 5.2 so that the diagrams in which the charged H + Higgs 
runs in the loop will be suppressed, to better explore the neutral Higgs sector properties. As we 
can see, the BR is highest when m h o is lowest, and rises with tan/3. The dominant diagrams in 
this case are the ones which have two h° or two H + scalars in the loop: The diagram with two 
h° dominates in the low tan (5 - low m h o region, while the diagram with two H + dominates in 
the high tan/3 - high m^o region, and is responsible for the rise of the BR with tan/3. Both of 
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(a) (b) 

Figure 5.1: (a) 3D plot of BR (t -> cH°) in the m H + - tan (3 plane in the T2HDM, and (b) the 
dominant diagram. We set m h o = 1000 GeV and m^o = 1200 GeV. The color scale represents 
the BR: the blue represents the lowest BR and red the highest. 




10 10TO u 

n 

(a) (b) 

Figure 5.2: (a) 3D plot of BR(t -> cH°) in the m h o - tan/3 plane in the T2HDM, and (b) the 
dominant diagram. We set m#+ = 1000 GeV and m A o = 1200 GeV. 
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Figure 5.3: 3D plot of BR (t — > off ) in the m^o — m/jo plane in the T2HDM. We set mjj+ = 
660 GeV and tan (3 = 28. 



parameters 


SM 


2HDM-II 


T2HDM 


m h o = 800, m A o = 1000, tan/? = 72, = 200 


6.03 x 10" 14 


4.25 x 10~ 5 


5.99 x 10~ 5 


m h o = 800, m A o = 1000, tan/3 = 72, = 380 


6.03 x 10~ 14 


1.79 x 10~ 6 


2.57 x 10" 6 


m h o = 200, m A o = 4000, tan/5 = 20, m H + = 1050 


6.03 x 10" 14 


5.15 x 10~ 8 


9.39 x 10- 5 


m h o = 200, m A o = 1000, tan/5 = 20, m H + = 1050 


6.03 x 10" 14 


3.34 x 10~ 12 


3.14 x 10~ 7 



Table 5.1: Comparison of the BR(t -> cH°) within the T2HDM, the 2HDM-II, and the SM. 
Masses are in units of GeV. 



these diagrams receive an enhancement from the 3-scalar vertex with large m^o , as mentioned 
above. 

The plots are similar, yet the BR are higher when the H + runs in the loop. That is a distinctive 
property of the T2HDM: the charged Higgs coupling H + bc can be enhanced by as much as 
Vtb/Vo, compared with any other 2HDM. 

In Fig. 5.3 we show the BR in the m A o — m^,o plane in the T2HDM, at the best fit of the 
charged sector parameters. The graph shows that the BR rises when m^o is highest and m^o 
lowest, since then the diagram with two m h o starts to dominate. The dip in the middle of the 
surface is due to cancellation in the H°H + H + vertex. 

In Figs. 5.4 and 5.5 we give 2D plots of the BR as a function of tan/3 and m H + respectively, 

with the same parameters as in Fig. 5.1. We can now see the dependence of the BR on tan/3 

and m H + more clearly: the BR rises with tan j3 but increases with smaller m H + . 

Finally, we wish to illustrate more clearly the difference between the T2HDM, the 2HDM-II, and 

the SM. For that purpose we give in table 5.1 the BR (t — > cH°) values within these 3 different 

models, for several points in the relevant parameter space. We recall that the 2HDM-II has a 

Yukawa potential similar to the MSSM, and has no tree-level FC interactions. 

The first two rows illustrate the impact of the charged sector, by setting a high m h o. The BR 



21 



CD 
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Figure 5.4: The BR(t — > cH°) as a function of tan/3 at various m H + in the T2HDM. We 
set m h o = 1000 GeV and m A o = 1200 GeV. "LHC thresh." stands for the limit of the LHC 
sensitivity at 100 fb^ 1 . 
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Figure 5.5: The BR{t — > cH°) as a function of m H + at various tan/3 in the T2HDM. We 
set m h o = 1000 GeV and m A o = 1200 GeV. "LHC thresh." stands for the limit of the LHC 
sensitivity at 100 fb^ 1 . 
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is a bit higher in the T2HDM than in the 2HDM-II. We note that the value m H + = 200 GeV 
is outside the la bounds. We can see that, for a high m^o, the BR(t — > cH°) in the T2HDM 
is not as enhanced as expected relative to the 2HDM-II. We expected that the diagram with 
two charged scalar and a b quark will be enhanced in the T2HDM, since the H + cb interaction 
is enhanced. We recall the Feynman rule for this interaction: 

-ji- — [tan (3V cb m b R + m c (- tan (5V ch + f * (tan (3 + cot (5) V tb ) L\ , (5.1) 
and the amplitude of this diagram: 



Mr = ^9 3 H h +H+h [m b C (a^B^L + B» + A» + *R) ~ m c C 12 (b» + B» + *L + A? A» + *r) + 
+m t (-C n + C 12 ) (A» + At*L + B^B^r)} «*, (5.2) 

where are the Passarino-Veltman scalar functions (see App. D). 

The term multiplied by the left projection operator is enhanced. In our notations it is denoted 
as A^ b + in (5.2). The leading term in the amplitude in the 2HDM-II, for tan/? > 10, is 
oc m t B^ b + B^ + * oc m t m 2 tan 2 (3V cb V t l (see Sec. 3). In the T2HDM this term does not receive a 
significant enhancement. However, the term oc m b A^ b + B^ + * , which in the 2HDM-II is a sub- 
leading term, is in the T2HDM oc £*m c m 2 tan 2 pV tb V t * b , and is of the same order of magnitude 
as the leading term of the 2HDM-II (together with the CV, loop functions). Thus we can 
summarize that what would have been a sub-leading term in the 2HDM-II, becomes in the 
T2HDM of the same order of magnitude of the leading term, and therefore the enhancement is 
not as significant as expected. 

In the last two rows of the table we see the impact of the neutral Higgs sector, by setting a 
high rriH+- The results are much higher in the T2HDM than in the 2HDM-II, which is to be 
expected, since the 2HDM-II does not have any tree-level FC interactions. However, the overall 
BR's in the small rrtyfi regime are also small due to the cancellation of the leading term in the 
h°it vertex. 

To complete the picture, we note that the SM value is only dependent upon the neutral Higgs 
mass. Setting m H o = 91 GeV, we get: BR SM (t -> cH°) = 6.03 x 10~ 14 , for m H o = 100 GeV 
we get BR SM (t -> cH°) = 4.63 x 10" 14 , and for m H o = 150 GeV we get BR SM (t -> cH°) = 
5.26 x 10" 15 . 



5.2 Results for the 1-loop Higgs rare decay H° —> tc 

We recall the values that we use in the following plots, which were also given above: m H a = 
300 GeV (chosen arbitrarily), m t = 172.5 GeV, m c = 1.24 GeV, m b = 4.20 GeV, m w = 
80.40 GeV, m z = 91.188 GeV, a EW (m z ) = 1/127.9, and a = 0. The f and f parameters 
are set to their best-fit value of Eq. (4.1): |f | = 0.8, ^ = 110°, \£\ = 0.21, and <p? = 250°. 
The total Higgs width is calculated (see Sec. 4.4) from the decays H° — > qq, H° — > VV, 
H° — > hihj and H° — > Vihj, as defined in App. E. 

In Fig. 5.6 we present the SM value for the BR(H° ^ic + ct). Our results agree with [4]. 
Next we turn to results in the T2HDM. In Fig. 5.7 we give a 3D plot of BR (H° ^ic + ct) 
in the m H + — tan/3 plane in the T2HDM. We see the same tendency as in the decay t — > cH°: 
The BR rises with tan/3 and rises with lower m H +. 

In Figs. 5.8,5.9 we give 2D plots of the BR as a function of tan/3 and m H + respectively, with 
the same parameters as in Fig. 5.7. We see now more clearly the behavior described above. 
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Figure 5.6: The SM value for the BR(H° — > tc + ct) as a function of the Higgs mass, for 
mb(Wb) = 4.2 GeV and for m&(mz) = 3 GeV [27]. The BR is not sensitive to m c . 




(a) (b) 

Figure 5.7: (a) 3D plot of BR (H° ^ic + ct) in the m H + - tan/3 plane in the T2HDM, and (b) 
the dominant diagram. We set m^o = 1000 GeV and m^o = 1000 GeV. 
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Figure 5.8: The BR (H° ^tc + ct) as a function of tan/? at different m H + in the T2HDM. We 
set m h o = 1000 GeV and m A o = 1000 GeV. 
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Figure 5.9: The BR (H° — >• tc + ct) as a function of mu+ at different tan (3 in the T2HDM. We 
set m h o = 1000 GeV and m A o = 1000 GeV. 
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(a) (b) 

Figure 5.10: (a) 3D plot of BR (H° -> ic + ct) in the m h o - tan/3 plane in the T2HDM, and 
(b) the dominant diagram. We set m H + = 1000 GeV and m A o = 1000 GeV. 



parameters 


SM 


2HDM-II 


T2HDM 


m h o = 800, m A o = 1000, tan/? = 72, m H + = 200 


1.23 x 10~ 13 


1.26 x 10" 4 


1.70 x 10" 4 


m h o = 800, m A o = 1000, tan/3 = 72, m H + = 380 


1.23 x 10~ 13 


3.09 x lO" 6 


4.45 x 10^ 6 


m h o = 200, m A o = 4000, tan/3 = 20, m H + = 1050 


1.23 x 10" 13 


8.69 x 10~ 8 


2.90 x 10~ 4 


m h o = 200, m A o = 1000, tan/3 = 20, m H + = 1050 


1.23 x 10" 13 


8.99 x 10~ 12 


9.11 x 10~ 7 



Table 5.2: Comparison of BR(H° -> tc+ct) between the T2HDM, the 2HDM-II, and the SM. 
Masses are in units of GeV. We set m H o = 300, a = /3, and other parameters to their best-fit 
value of (4.1). 



In Fig. 5.10 we give a 3D plot of BR (H° ->tc + St) in the m h o - tan/3 plane in the T2HDM, 
and in Fig. 5.11 we give a 2D plot of the BR as a function of tan/3 with the same parameters 
as Fig. 5.10 at several values of m ft o. We can see that the BR rises with lower m h o, but has a 
weak dependence on tan/3. 

In table 5.2 we give the BR (H° — > ic + ct) values in the different models, for a few points in 
the parameter space. As can be seen, the behavior is similar to the t — > cH° process, although 
generally the BR values are higher. 
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Figure 5.11: The BR (H° -> tc + ct) as a function of tan (3 at different m h o in the T2HDM. We 
set m H + = 1000 GeV and m A o = 1000 GeV. 
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Chapter 6 
Summary 



The T2HDM is a distinct type of a 2HDM in which the top quark receives a special status. In 
this model, the top is coupled to the second Higgs doublet, while all other quarks are coupled 
to the first Higgs doublet. Assuming that the second Higgs VEV (1*2) is much larger than the 
first Higgs VEV (i>i), the top quark receives a much larger mass than all other quarks in a 
natural manner. Therefore the working assumption of the T2HDM is that tan/5 = v^/vi 
1. In addition, these Yukawa couplings generate potentially enhanced flavor-changing (FC) 
interactions, both in the charged and the neutral sectors. These interactions can greatly enhance 
FC decays such as t — > ch and h — ■> ic. 

The Yukawa sector of the model was explicitly (and independently) derived, as well as scalar 
self interactions. For example, the H + bc vertex is enhanced by a factor of V^/V^, compared 
to the corresponding 2HDM-II vertex. This enhancement motivated the present work, since it 
is expected to influence the 1-loop t — > cH° and H° — > ic decays via diagrams involving H + 
scalars and b quarks inside the loop. 

In order to separate the 1-loop decays from the tree-level decays, we chose a = j3. This choice 
eliminates the t — ► cH° and H° — > ic tree-level decays, so that these decays proceed at 1-loop. 
On the other hand, the decays t — > ch° and h° — > ic occur at the tree-level for a = f3. For these 
tree-level decays we gave explicit formulae. The h°ic neutral FC interaction can also enhance 
the 1-loop BR{t — > cH°) and BR(H° — > ic), via diagrams involving h° scalars and t quarks 
inside the loop. 

The parameter space of the T2HDM was explored for the resulting 1-loop BR (t — > cH°) and 
BR (H° — > ic). We focused on those regions of the parameter space in which these BR's can be 
much higher than in the SM and 2HDM-I,II. We found the dynamics of the two processes to be 
similar, which was expected since their amplitudes and rates are related by crossing symmetry. 
The 1-loop BR(t -> cH°) can reach ~ 10" 4 in the T2HDM. This is above the LHC detection 
threshold of 5 x 1(T 5 , and above the SM, and 2HDM-I,II predictions. The 1-loop BR (H° -> ic) 
can reach above ~ 10~ 4 in the T2HDM, higher than the SM and 2HDM-I,II predictions. 
We conclude that if such decays are indeed identified at the LHC, then the dynamics of the 
T2HDM type will be especially motivated. 
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Appendix A 



Higgs potential in two Higgs doublet 
models 



In the following we introduce the Higgs potential of a general 2HDM, which corresponds to the 

T2HDM as well as to 2HDM's of types I, II and III. 

We assume a CP conserving Higgs potential of the form [5]: 



+A 4 



(<^$ 2 ) - ($I$ 2 ) ($2$l) + ^5 $I$2 - UlV 2 



+ 

(A.l) 



already introduced in Sec. 2, and: 



$1,2 = 



<^2 

t)l,2+*? r 2 +i' 



? 1,2 



(A.2) 



This potential has five couplings A» plus two VEV's i>i and v 2 , seven degrees of freedom in total. 
These will be later expressed in terms of 5 masses of the physical scalars, plus two angles. We 
will then extract the Feynman rules of the 3-scalar interactions, expressing them in terms of 
the physical masses and angles. 

We assume a potential which conserves CP. The absence of CP violation implies that the CP- 
even h°,H° and the CP-odd A°,G° mix separately, as we will later see. CP violation in the 
scalar potential would mix CP-even and odd scalars (H°, h°, A ), as discussed in [5]. 
The wearying part is deriving the Feynman rules in terms of the masses and angles instead of 
Aj ,vi and t> 2 . The explicit derivation is straightforward and we will not follow it completely 
here. We will, however, introduce the important formulae, following the notation of [5]. 
The fields can always be redefined so that their VEV's are real, without affecting the potential. 
The VEV's as defined can be easily seen to minimize the potential. 

The mass terms of the neutral real (CP-even) scalars can be combined into a symmetric bilinear 
mass term: 



1 / Av\ (Ai + A 3 ) + w|A 5 (4A 3 + A 5 ) v x v 2 

2 1 (4A 3 + A 5 ) v x v 2 Av 2 2 (A 2 + A 3 ) + v*\ 5 



The mass-squared matrix, M, can be diagonalized using the rotation matrix 
with an angle a such that: 



Or 



(A.3) 



cos a -sintt 
sin a cos a 
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sin 2a = 



2M 



12 



cos 2a = 



M n - M 22 



y/(M n -M 22 ) 2 + 4M 1 2 2 ^(Mn -M 22 ) 2 + 4M 1 2 2 

which defines the CP-even neutral scalars in the mass basis: 



(A.4) 



H° = $f cos a + $f sin a, 
h° = -$° r sin a + $ 2 r cos a, 



with masses: 



h° = o ( M H + M 22) ± o V ( M H - M 22) + 4M ] 



[ 12 



The CP-odd neutral scalars have a simpler mass 2 matrix: A5 
ized using the rotation matrix with the angle f3 such that: 

tan/5 = — , 

Vi 



-ViV 2 



-ViV 2 v{ 



(A.5) 



(A.6) 
diagonal- 

(A.7) 



which defines the CP-odd neutral scalars and unphysical Goldstone boson (G°), in the mass 
basis: 



with mass: m 2 A0 = A 5 w 2 /2 (recall that: v = \J v \ + 



G° = $f cos/3 + $£ sin/5, 
A° = -$fsin/3 + $™cos/3, 



(A.8) 



We can see that the CP-even h°, H° and the CP-odd A , G° mix separately, as promised. 



The charged scalars have a similar mass 2 matrix: A4 



2 \ _ 



viv 2 



-viv 2 

v 2 



, diagonalized with 



the same angle f3, defining the charged physical scalars and unphysical Goldstone bosons (G^), 
in the mass basis: 



G ± = $f cos[3 + $f sin/3, 
if± = -$fsin/3 + $£ cos (3, 



(A.9) 



with mass: rn 2 H± = X4V 2 /2. 

As we stated above, we want to express the triple-scalar couplings in terms of masses and angles 
instead of \i,Vj. We will first connect between and \,Vj\ 



( M n ^ 




' Av\ 





Av 2 


M 22 


\- 





Av\ 


Av 2 


V M 12 J 










Aviv 2 




(A.10) 



We can invert the equation, and insert A5 = Hr-, to write: 
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Ai + A 3 = - —3- - tan /3^- 
4 y f j v 2 J 

A 2 + A 3 = - —j- - cor /3^- 

2 4 \t>it>2 f / 

The neutral CP-even couplings can be collected from (A.l): 



£CP-even = 4Wl ($*) 3 (A X + A3) + 4«a ($?") 3 (A 2 + A3) + 

+4 Wl $r ($r) 2 (a 3 + ^A 5 ) + An*? «) 2 



A 3 + ^A 5 



3 / M n . 2 ^ m ^o 



f 1 (Hq cos a — h sin a) — 5 tan /3 — ^- + 

(iW 771 2 

— 5 cot 2 /3 — 4" ) + 
v£ v 2 1 



Js/L tHj^ 

+vi (H cos a — h sin a) (Hq sin a + h cos a) 2 ( 1 4~ ) + 

1 ' 



+t>2 (i?o sin a + ho cos a) (Hq cos a — h sin a) 2 ( 1 ^ ) . (A. 12) 



For example, focusing on the h°h°H° term, and using v\ = v cos/5, v 2 = v sin/5, we get: 



If „ ■ 2 I . 2 o 2 \ , o • o • 2/ ^22 ia 2 



£h°h°H° = _ { 3 cos /3 cos a sin a — — — tan (3m A o + 3 sin (3 sin ct cos a — t, cot (3m A0 + 

v { \cos 2 (3 J \sin (3 J 

+ [cos (a + /3) cos 2a — sin (a + (3) sin ct cos a] ( — — — - + m 2 A0 ) 1 . (A. 13) 

L v y v y J \cos/3sm/3 /J v y 

In order to arrive to the final form, more algebraic work is needed. Using (A. 4), from which 
follows also: 2M i2 cos 2a = (Mu — M22) sin 2a, and (A. 6), from which follows also Mn + M 22 = 
m h° + m H°i an d other trigonometric identities, we arrive at: 



£h°h»H° = ^ . n J [sin 2a (2m 2 + m 2 H0 ) - m\ (3 sin 2a - sin 2(3)] , (A. 14) 

Zmy/ sin Zp 

which agrees with [15], and which is expressed in terms of the physical masses and mixing 
angles. 

All other 3-scalar interactions are derived similarly. In App. (B) we listed all the vertices 
relevant to the present work. 

We note that we did not find the h°H°H° coupling for the higgs potential of (A.l) in the 
literature and it was therefore derived by us. 
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Appendix B 

Feynman rules for two Higgs doublet 
models 

In this section we list the Feynman rules for the 2HDM's that were used in this work, in the 
t'Hooft Feynman gauge. The Feynman rules are presented in Fig. B.l. 

In table B.l we give the Yukawa Feynman rules of the T2HDM and of the 2HDM type II. The 
other Feynman rules are common to all 2HDM's. 

In table B.2 we give the vector- vector-scalar couplings, common to all 2HDM's, from [5]. 

In table B.3 we give the vector-scalar-scalar couplings, we define the vertices as in [5], where 

the second particle is outgoing. 

The vertices Z°G°H°, Z°G°h° do not participate in the calculations since the corresponding 
Yukawa vertex qqG° does not generate FC interactions. 

Now we turn to the 3-scalar interactions. We were not able to find all vertices in the literature, 
and therefore we derived one vertex (h H°H ), while the rest of the scalar self interactions can 
be found in [15]. We give in table B.4 the complete list of the 3-scalar interactions that were 
used in this work. 
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Figure B.l: Feynman rules. 



T2HDM 



2HDM-II [5] 



H°UjUi 



2mw 



u cos p 



+ 



sin p cos 



R+(h.c.)L 



gM u sina 
2miy sin/3 



h° 



UU 



2m, 14/ 



" cos p 



sin P cos /3 



R+(h.c.)L 



gM u cos a 
2mw sin/3 



' 2mw 



A 



MM 



2miy 



-M u tan /3 + E (tan /3 + cot (3)) R + (h.c.) L 



- cot P(R-L) 



G° 



UU 



(R — L) 



gMd cos a 
2mw cos P 
gMd sin a 
2mw cos 



>-(R-L) 



gMd cos a 
2mw cos /3 



gMd sin ct 
2m-i4/ cos P 



ij£ tan P(R-L) 



A°dd 



~g~XTd 



G°dd 



' 2m 



ill 



(i2-L) 



(R-L) 



H+ud 



[tan f3V CKM M d R 
M u tan/5 + E (tan /3 + cot V C kmL] 



y/2m\\, 

+ (- 



+ cot/3M M \/ CifM ^] 



G + ud 



\f2mw 



(M u Vc 



CRM 



L — Vt 



CRM 



M d R) 



V2mw 



(M U V C 



CRM 



L — V, 



CRM 



M d R) 



Table B.l: Feynman rules for Yukawa interactions in the T2HDM and in the 2HDM-II. 
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W + W-H° 


igmw cos (P — a) g^ u 


W+W-h" 


igmw sin (P — a) g^ v 


Z°Z°H° 


t9n l z cos (p a) g^ 

cos tiw ' 


Z°Z°h° 


19 ™ Z sin ((3 a) g^ 

cos Vw y ' 



Table B.2: Feynman rules for vector- vector-scalar interactions [5]. 



y a ^ 


= ig$* (P b + P c f 


W + H+H° 


if sin (f3- a) {P H+ +P HO f 


W + H + h° 


-if cos {P-a) (P H ++Ph°y 


W+G+H" 


-zf cos (P-a) (P G+ +P H »y 


W+G+h" 


-if sin (P- a) (P G+ +P h oY 


Z°A°H° 


-"gS? (Pa* + Ph*)" 


Z°A°h° 





Table B.3: Feynman rules for vector-scalar-scalar interactions [5]. 



H+H-H 


9_ 


(m^ + — m 2 A0 + \m 2 H0 ) cos (/3 — a) + (m^ — m Ht>) cot 2/3 sin (/3 — a) J 


H + H~h° 


g 

m w 


(m 2 H+ - m^ + \m 2 h0 ) sin (/3 - a) + [m\ - m^ ) cot 2/? cos (/3 - a) J 


h°h°H° 


LJL20 [( 2m lo + m Ho) sm 2 « m % (3 sin 2a sin 2/3)] 


h°H°H° 


LjLw [i 2m H° + m l°) sm 2 « m Ao (3 sm 2a + sm 2/3)] 




9 

2m, 14/ 


m^o cos (P — a) + 2 (m^ — co ^ 2/^ sin (/3 — a) 


A°A°h° 


2miv 


m 2 hQ sin (/3 — a) + 2 (m^ — m^ ) cot 2/3 cos (/3 — a) J 


H + G~H° 




H+G-h° 


l i ( m H+ m2 )cos(/3 a) 


G + G~H {) 


- l 2 ra >ff» COS (^ °0 


G + G~h° 





Table B.4: Feynman rules for triple-scalar interactions [15, 5]. 
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Appendix C 

1-loop diagrams calculation 



In this appendix we give the 1-loop calculation of the 10 diagrams shown in Fig. 4.1. The 
calculation was done in the t'Hooft Feynman gauge. 

In the t'Hooft Feynman gauge the vector bosons propagators reduce to their simplest form: 
A = ig fJ,v [p 2 — m 2 + ie] \ and the Goldstone bosons mass is set equal to the respective gauge 
bosons: m G + = m w +, m G o = m z . The t'Hooft Feynman gauge was chosen because the 
calculation of each diagram is simpler. 
definitions: 

M n - the amplitude corresponding to diagram n 
h - the external neutral scalar 

i - (— t) when used as index, the incoming fermion - the top 

j - (= c) when used as index, the outgoing fermion - the charm 

a,P - when used as indices, internal bosons (vectors or scalars) in the loop 

l,k,q - when used as indices, internal fermions 

L,R the Left, Right projection operators 

Uj - (= u(Pj) ) the outgoing spinor of the charm 

Ui - (= u(Pi) ) the incoming spinor of the top 

B , E>i, Co, Cij - the n-point integral functions, defined in App. D 

A s ab , B 5 ab - the left, right -handed parts of the fermion-fermion-scalar vertex 

a abi b S a b ~ the left, right -handed parts of the fermion-fermion- vector vertex, for both charged and 

neutral gauge bosons 

g3^vhh,wh ^ e ver f- ex f 3_ sca i arS) vector-scalar-scalar, vector-vector-scalar, respectively 
- the metric, g^ u = diag(l, —1, —1, —1) 
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Ill ■ — 1 

Mi = Y^mj^m 2 t mimkB ° ( S §^ S * L + A i! B ^ R ) + 



- m l m l B 1 (B?*A? k A? k *L + A^ B? k B? k * R) + mi m k B (B%B«A%L + A%A« k B? k *R) + 
-m 2 B, (B%B? k B? k *L + A%A%A$R)] u h (C.l) 

where B = B (m\, m 2 a , m 2 ) . 

ffl ■ 1 

M2 = Y^^p^ [mWkBo (A%B?:b£L + B^A^R) + 

+ m k m 3 B (Bj k Af k *B^*L + A%B% A* R) - m^B, (B^B^B^L + A%A%A%R) + 
-m)B x (A%A™BtL + B%B%A%R)] u t , (C.2) 

where B = B (m k , m 2 a , m 2 ) . 

Ms = Y^^fT^ [^m k B (B^bf^L + A%a? k b%R) + 

+ 2m l m i B 1 (B^b^L + Afta&agR) + Am t m k B (B%a? k b%L + Ap«a%R) + 
+2m 2 B 1 (B?*a? k a%L + A%bfcb%R)] u h (C.3) 

where B = B (m\, m 2 a , m 2 ) . 

Ill ■ ! 

M, = [±rn imk B (b%a%B^L + a%b%A%R) + 

+ 4m k m,B (a%b%B%L + b%a%A%R) + 2m j m l B 1 (a^B^L + fyfftApR) + 
+2m 2 B 1 (b%b%B%L + a^A^R)} u t , (C.4) 

where B = B (rn\, m 2 a , m 2 ) . 

M " = tB ^ L + { [ do + m * Cu + ^ - ^ ° 12 \ ( A ^< B % L + B n B ^R) + 

- m q m t Cn (A« q B h k * q Af k * L + B^B^R) + m q m 3 C l2 (B^At^L + A^AgR) + 

+ mi m 3 (C 12 - C n ) (B« q B%A? k *L + A^A^B^R) + m q m k C (A« q B^B? k * L + B^A^R) + 

- m mk (C n + Co) (AJ q A k *Af k *L + B^ q B kq B^ k *R) + 

+ mj m k (C 12 + C ) (B? q B£* q B? k *L + A%A h k * q A^R) } u h (C.5) 
where C = C (m\,m 2 a ,m 2 ,m 2 ,m 2 ,mfy . 
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Mt 



16n 2 



4C + 2 (m 2 - mj + m 2 ) C n + 2 (-m? + m 2 + m*) C 12 ] {b%B%a%L + a%A h k * q b% R) + 

+ 2m q m l C 11 {b%A h k * q b%L + a%B h k * q a%R) - 2m q m 3 C 12 {a%B h k * q a%L + b« q A h k * q b% R) + 

+ Am q m k C (b« q A h : q atL + a%B h k * q b^*R) + 2 mi m k (C n + C ) {b%B h k * q b% L + a« q A h k * q a%R) + 

-2m jmk (C 12 + C ) {a^A^L + b%B h k * q b%R) } u h (C.6) 



where C = C (m 2 k , m 2 a , m 2 q , m 2 , m 2 , m 2 h ) 



M 7 



-iu 



3 „3h 



16tt 2 



9a/3h 



m k C [A%B%L + BP jk AgR) - m,C 12 [BP jk B£L + A%A%R + 



+m t (-c n + c 12 ) (a^a-l + fi&i^i?)' 



(C.7) 



where C = C (ml, m^, m 2 ^, m 2 , m 2 ,, m 2 ^ . 



16tt 2 



9a{3h 



4m fe C (b%a%L + o^fl) + 2m, (C n - C 12 ) (b%b%L + a^fl) 



+2m,C 12 (aJa^L + ^R 
where C = C (m 2 k , m 2 , m 2 p , m 2 , m 2 h , m 2 ) . 



U.; 



+ 



(C.8) 



Mo 



-vhh 



lQ n 2 9 Pah 



Co + 2m 2 C n + mjC 12 - 2m 2 C 12 (b^L + aj^fl + 



mi m 3 (C 12 + C n ) (aJ^L + ^^i?) + m,m fc (C - C 12 ) (a^L + ft^i?) + 



+m J m fe (C 12 - C n - 2C ) (&^ fe *L + a%B^R) 
where C = C (m 2 k , m 2 a , m 2 p , m 2 , m|, m 2 ) . 



(C.9) 



Mi 



10 



16tt 2 



3af3h 



m 2 (C 12 - C n ) - 2mJC u - 2m* (C 12 - C n )) (^aJfL + fif^i?) + 
+ m z m 3 (2C n - C 12 ) (b^L + A%agR) + m 3 m k (C 12 + 2C ) (s&a£L + A%b£R) 4 

k, (CIO) 



+m,m fc (C n - C 12 - C ) (Aj^L + B&a£i2 
where C = C (m 2 k , m 2 a , m 2 p , m 2 , m|, m 2 ) . 



37 



Appendix D 

Definition of the n-point integral functions 



We present here the definitions for 1-loop scalar, vector and tensor integrals: 

Bo;BpW , f4p ^/g ;^ (D ,) 



C ; C M ; C^; C (m 2 , mj, m\ pj, pf) 



d^k 1; k^~, k^ u ] k 2 



in 2 [k 2 - m\\ [(k + p^ 2 - m 2 ] [(k + p 1 + p 2 f - ml] ' 

(D-2) 



B M =P IJ ,B l , 

CW = PimPi^ C 2i + P2fiP2vC 22 + {piP2}^C 23 + g^C 24 , (D.3) 



where {afr}^ = a M &„ + a„6 M . 
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Appendix E 

Higgs width calculation 



In this section we give the formulae that we used in calculating the Higgs width. The main 
contributions to the total width Y tot are: 



^tot 



Y^h^qq _|_ pfo— ►W _|_ Y^~*HiHj _|_ p/i— *VH 



Only leading order values were used. The decay products are all taken to be on-shell, and their 
secondary decay products are not taken into account. Each contribution of h — > x + y was 
calculated above the threshold: m h > m x + m y . The vertices are defined here as in Fig. B.l. 
The leading order value of h — > qq is [5]: 



NcA\ m ( imp i 



r^H-g^i--^ , (E.2) 

where A hqq = — J^^f ; ^m^craf is the quarks coupling to H°; h°, respectively, and N c = 3 is 
the color factor, as mentioned above. We also give explicitly the width for the process h° — > bb 
in the T2HDM for a = (3: 



3g 2 ml 2 f Am 2 



The leading order value of h — > W + W~ is [5]: 

r (h - w + w~) = 9 ^w< {l _ x) \( l _ xJr l x 2\ (E4) 

V ; 64vrm^ K ' \ 4 J ' V ' 

where gnww =9 r nw cos (P ~ a )', gmw sm (P ~ a ) is the W + W~ coupling to H°; h°, respectively, 

ml 

The leading order value of h — > Z°Z° is [5]: 

r(^^-&rf^(i-.)»(i-. + ^), (E,) 

where = cose w cos — a )' cos^7 s * n — a ) * s t ne coupling to H°; h°, respectively, 

, 4m| 

and x = — f. 

Since the couplings W + W~h° and Z°Z°h° are both oc sin (P — a) (as we have shown in table 
B.2), then by choosing a — P, the widths of the processes h° — > W + W~ and h° — > are 
both reduced to zero. 
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Figure E.l: The total width of the SM Higgs: leading order approximation compared to cor- 
rected width of [16]. 



The leading order value of h — > HiHj (where Hi and Hj are any two scalars) is: 

^^Hf^4't# (e - 6) 

where ghUiH is the triple Higgs coupling hHiHj. 

The leading order value of h — > VH (where VH = W + H~ or Z° + neutral scalar) is: 

V f h , y H \ _ 9vH h ™ 2 v x i (i m v m H\ X (l m h m H\ ( E7 ) 

where guvH is the vector-scalar-scalar coupling VHh. This contribution can be important, as 
stated in Sec. 4.4. 

In order to demonstrate the error introduced by including only leading terms in the Higgs width 
calculation, we give in Fig. E.l the total SM Higgs width as calculated in this work compared 
to the width of [16] which includes higher-order corrections. 

As can be seen, below the WW threshold (at about 160 GeV) the values are different. In this 
mass range bb decay dominates, and corrections have a large impact. However, this mass range 
is also below the h — > tc threshold. On the other hand, above the WW threshold the values are 
very similar, and corrections have a small impact, and so the use of leading order approximation 
can be justified. 
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Appendix F 



Cancellation of divergences in the 1-loop 
amplitude 



Some of the 1-loop Feynman diagrams have a divergent part. These divergent parts cancel, 
since this is a leading-order calculation and, therefore, there is no renormalization. 
This cancellation is also important as a means of checking the self consistency of the calculation. 
We define e — 4 — d, where d — > 4 is the number of dimensions. As e — > 0, some n-point 
integrals will have a term proportional to -. These are summarized below: 



Bn 



€ 

e 



C*24 ~ 
C ~ - 



1 1 

~2? 

-2 1 -. 

e 



(F.l) 



The parts of the Feynman diagrams proportional to ^ were collected below. Only the parts of 
the diagrams with a left projection operator are given, while the right-handed parts are subject 
to a similar cancellation. The Feynman diagrams were shown in Fig. 4.1. 



M. 



ool=j 
1L 



1 



m: 



ool=i 
2L 



mf 


-mf 




-1 


m) 


-mf 




1 


mf 


-mf 




1 






m) 


-mf 



—2m/ 
- rrii 
-2m k 

-rrij 
-8m k 

— 2mi 
-8m k 

— 2m i 



rniB»°*A? k B? k * + m^B^) - 
mB^Af^ + mtB^B-B?:)' , 
mA^B^Bf* + mjBf^B^) - 

rniB%B?:B»°* + m,A%^Bf *)' 
miB»°*b? k a% + m^f *«*) - 
mtB^bX + m i B»°*a? k a° k *)~ , 

m l b jk a lk B il + m i a jk b lk B il ) ~ 

mia jk a lk B u + mjb jk b lk B u J 
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M oo q = k = 2A'> q A»°*B%, 

M? L = 0, 
= 0, 

M? 0L = 2g$> H0 A%a£. (F.2) 

After inserting the Feynman rules of the T2HDM, we were able to show that the terms pro- 
portional to - cancel as shown below: 



Mj 30 


+ M 2 °° 


+ M 5 °° 


M3 00 


+ M 4 °° 


= 0, 


M 6 °° 


= 0, 




M 7 °° 


= 0, 




M 8 °° 


= 0, 




M 9 °° 


= 0, 




MS 


= 0. 





(F.3) 



This cancellation was also verified numerically in the FORTRAN code. 
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